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Abstract 

Time-dependent Schrodinger's equation is integrated for a one-dimensional strongly-correlated electron system 
driven by large electric fields. For larger electric fields, many-body Landau-Zener tunneling takes place at anti- 
crossings of the many-body energy levels. The nonlinear I-V characteristics as well as the time dependence of the 
energy expectation value are obtained. The energy of the Mott insulator in electric fields shows a saturation, which 
suggests a dynamical localization in energy space of many-body wave functions. 
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Introduction Strongly-correlated electron systems 
should be an interesting test-bench for nonequilibrium 
phenomena. Specifically, we can ask how the nonlinear 
transport properties of Mott insulators are distinct 
from those of band insulators, ft has been experi- 
mentally shown that a quasi-one-dimensional cuprate 
exhibits a dielectric breakdown that remains when ex- 
trapolated to zero temperature[l]. This suggests that 
the origin of the breakdown is quantum mechanical. 
Indeed, a many-body version of the Landau-Zener 
transition for correlated electron systems has been 
proposed by the present authors[2]. In the present 
paper we make a further analysis of the many-body 
Landau-Zener breakdown of the one-dimensional Mott 
insulator in strong electric fields to elucidate how its 
nature differs from those of the conventional Zener 
breakdown of band insulators. 

I-V characteristics We study the time evolution of 
wave functions for a one-dimensional Mott insulator in 
a static electric field F turned on at t — 0. The electric 
field is here induced by a time-dependent Aharonov- 
Bohm(AB) flux (p(t) = $(*)/$o = FLt/h pierced 
through the system in a periodic boundary condition, 



so the evolution is governed by the time-dependent 
Schrodinger equation with a time-dependent Hamilto- 
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When the system is half-filled, the groundstate is a 
band insulator for the level offset A > with the Hub- 
bard (7 = 0, and becomes a Mott insulating when U > 
with A = 0. 

We first look at the I-V characteristics (Fig. I) of the 
Mott insulating phase (U > 0, A = 0), where the cur- 
rent is defined as a time-averaged (J) = ^ f Q {J(t))dt. 
The current is suppressed until the field strength ex- 
ceeds a threshold, F(U), whose value depends on U. 
The threshold can be estimated from the Landau- 
Zener formula, and we predict an F(U) oc [AE(U)] 2 
dependence on the charge gap AE(U) as confirmed 
numerically [2]. Above the threshold, Landau-Zener 
tunneling from the ground state to excited states is 
activated. 
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Fig. 1. I-V characteristics for the one-dimensional Hubbard 
model for an N — 6 site ring. 

What distinguishes the breakdown of Mott insula- 
tors from those for band insulators is how the system 
becomes out of equilibrium in electric fields above the 
threshold. In order to quantify this, we plot the long- 
time behavior of the energy expectation value (H(t)) 
for a Mott insulator as compared with that for a band 
inslator(Fig.2). 

Small F regime: When the Landau-Zener tunnel- 
ing is ineffective, the wave function approximately 
evolves adiabatically along the ground state of the 
time-dependent Hamiltonian H(<f>(t)). The ground- 
state energy shows an AB-oscillation for a finite 
system. 

Large F regime: (H(t)) of band insulators diverges 
in the long-time limit, which is not apparent in 
Fig. 2(a), but the long-period oscillation seen in the 
figure has a periodicity A(/> ~ N, so the amplitude 
as well as the periodicity diverge for N — > oo. This 
is in sharp contrast with the saturated behavior of 
(H(t)) of Mott insulators. 
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Fig. 2. Time evolution of the energy (H(t)) of a band insula- 
tor(a) and a Mott insulator(b). 

The divergence of (H(t)) in the band inslator can 
be readily understood. The Landau-Zener tunneling 
produces free electron-hole pairs, which are accelerated 
by static electric fields. Since there is no dissipation for 
the present clean and isolated system nor electrode to 
absorb the energy, the energy increases indefinitely. 

On the other hand, (H(t)) in the Mott insulator is 
seen to saturates to a value that depends on the elec- 
tric field, which is quite an interesting and nontrivial 
behavior. For noninteracting but disordered systems, 
Gefen and Thouless explained such a saturation in 



terms of a dynamical localization, i.e., Anderson's lo- 
calization in energy space[3]. We can extend the notion 
of the dynamical localization to a (clean) many-body 
system. [4] This takes place on many-body adiabatic 
spectra typically depicted in Fig.3. We can see that, 
while a band insulator has a single level anticrossing 
separated by a gap between the ground state and the 
state excited with one electron-hole pair, a Mott insu- 
lator has many (actually a macroscopic number of) an- 
ticrossings among the excited levels. At each level an- 
ticrossing the Landau-Zener tunneling takes place and 
the probability amplitude bifurcates. 

An important point to note is that this branching 
occurs in a quantum mechanical manner, where con- 
tributions, with their phase factors, from various paths 
(in the energy-level vs time space) interfere with each 
other, resulting in the dynamical localization. The am- 
plitude of the distribution resides near the ground state 
when the electric field is weak, whereas a component 
extending into excitated states arises when F becomes 
larger, which is governed by the above phase factor. 
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Fig. 3. The adiabatic level scheme for a band insulator(a) and 
a Mott insulator (b). The level anticrossings arc highlighted. 

Conclusion 

So, while the disorder gives rise to a saturated be- 
havior in the current-time characteristics in one-body 
systems, this is replaced by the quantum mechanical 
tunneling across many level repulsions caused by the 
electron-electron interaction in many-body systems. 
We have identified this as the reason why the strong 
field behaviors are quite different between band insu- 
lators and Mott insulators. 
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